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Abstract. The magnitude- redshift relation is one of the 
tools for a direct observational approach to cosmology. 
The discovery of high redshift Type la supernovae (SNIa) 
and their use as "standard candles" has resurrected inter- 
est in this approach. Recently collected data have been 
used to address the problem of measuring the cosmologi- 
cal parameters of the universe. Analysed in the framework 
of homogeneous models, they have yielded, as a primary 
result, a strictly positive cosmological constant. However, 
a straight reading of the published measurements, con- 
ducted with no a priori idea of which model would best 
describe our universe at least up to redshifts z ~ 1, does 
not exclude the possibility of ruling out the Cosmological 
Principle - and cosmological constant - hypotheses. It is 
therefore shown here how the large scale homogeneity of 
this part of the universe can be tested on our past light 
cone, using the magnitude-redshift relation, provided suf- 
ficiently accurate data from sources at redshifts approach- 
ing z — 1 would be available. An example of an inhomoge- 
neous model with zero cosmological constant reproducing 
the current observations is given. The presently published 
SNIa data can thus be interpreted as implying either a 
strictly positive cosmological constant in a homogeneous 
universe or large scale inhomogeneity with no constraint 
on A. An increase in the number and measurement ac- 
curacy of the candidate "standard candles" at very high 
redshift is therefore urgently needed, for progress in both 
fundamental issues of the Cosmological Principle and of 
the cosmological constant. 



dSandage 196Tj IKristian fc Sachs 19661 lEllis et al. 1985)) . 
First presented as an ideal program which could be 
realised "in principle" , it has of late become a new 
field in cosmology, due to recent tremendous progress in 
instrumental technology. 

The magnitude-redshift relation is one of the tools 
for a direct observational approach. The discovery 
of high redshift (z ~ 0.5 — 1) Type la supernovae 
(SNIa), and of their potential use as "standard can- 
dles" dPhillips 1993) |Riess, Press fc Kirshner 1995) 
Hamuy et al. 1995) H996al Il996bl Perlmutter et al. 
19971 |Tripp 1998| Uha et al. 1999)) has resurrected in- 
terest in this approach. Recently collected data, from 
ongoing systematic searches, have been used to address 
the problem of measuring the cosmological parameters of 
the universe. Analysed in the framework of homogeneous 
models, they have yielded, as a primary result, a strictly 
positive cosmological constant, many orders of magnitude 
smaller than the energy of the vacuum expected in 
standard particle physics models, as proposed by Riess et 
al. (|1998[) . hereafter refered to as R98, and by Perlmutter 
et al. (|T9"99"1) . hereafter refered to as P99. 

A non zero cosmological constant was long ago pro- 
posed as a possible interpretation of magnitude-redshift 
observations of cluster galaxies and quasars by Solhcim 
( 1966J) . But the accuracy of the data was not then 
sufficient to probe this hypothesis. 
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1. Introduction 

The use of astronomical observations to di- 
rectly determine the space-time geometry of our 



universe 



is 



a long time proposed instrument 



Correspondence to: celerier@obspm.fr 



Anything that contributes to the energy density of the 
vacuum acts just like a cosmological constant. Analysed 
in the framework of Friedmann models of the universe, 
this can be viewed as an argument in favor of a zero 
cosmological constant, since some unknown symmetry 
of particle physics theory could presumably cancel the 
vacuum energy density (see e.g., |Weinberg 1989). If the 
value of the cosmological constant was confirmed to be 
in the range favored by the SNIa announced results, it 
would be necessary to explain how it is so small, yet non 
zero. As such a result would have a revolutionary impact 
on our understanding of the fundamental laws of physics, 
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it is of tremendous importance to check the different 
available interpretations. 

Provided every other source of potential bias or 
systematic uncertainties has been correctly taken into 
account, one of the most vulnerable points of the data 
collecting procedure is the difference in the absolute 
magnitude of the supernovae due to an evolution of 
their progenitors. This point is discussed at length 
in R98 and P99 who conclude that their results can- 
not be markedly affected by this effect. Other authors 
flDrell, Loredo fc Wasserman~ 1999, Domi ngez et al. 1999D 
claim contrary evidence that high and low redshift su- 
pernovae observed so far, substantially differ from one 
another, and thus contest their use as "standard candles" . 
A progressive dimming of the SNIa by intergalactic dust 
has also been proposed as a possible systematic effect that 
could mimic the behaviour of cosmic acceleration, and 
thus allow cosmologies without a cosmological constant 
(Aguirrc 1999). This important problem is not treated in 
the present work, where the working assumption is made 
that the SNIa data actually measure the magnitude- 
redshift relation up to the precisions claimed in R98 and 
P99. 

The theoretical interpretation is another matter which 
has to be dealt with carefully to avoid a priori assertions 
which would lead to incorrect results. 

P98 discuss the possibility that the magnitude- 
redshift relation they find from the analysis of their 
data is due, not to a cosmological constant, but to 
an evolving field of unknown nature that contributes 
to the total energy density of the universe (see e.g., 
ISteinhardt 19961 |Caldwell, Dave fc Steinhardt 1998[ 
IGarnavich et al. 1998[) . This entity, which can have an 
equation-of-state ratio different from that of the cos- 
mological constant, would lead to a different expansion 
history. 

It is proposed here to focus on a less exotic alternative, 
namely the possibility of large scale inhomogeneity of 
the part of the universe which can be probed with SNIa 
measurements . 

Large scale spatial homogeneity of the universe is a 
belief most commonly shared by current cosmologists. It 
procceeds from a hypothesis brought to the status of Cos- 
mological Principle by Einstein (|1917p . Its justification is 
based on two arguments: 



1. The isotropy, or quasi-isotropy, of the temperature 
of the Cosmic Microwave Background Radiation 
(CMBR) around us. 



2. The Copernican assumption that, as our location 
must not be special, the isotropy we observe must be 
observed identically from any point of the universe. 



Since a matter distribution which is seen isotropic 
from everywhere implies homogeneity, the Cosmological 
Principle follows. However, of the two above arguments, 
only the first is observation grounded. The second, purely 
philosophical, has never been verified, and cannot thus 
constitute acceptable evidence. 

As will be shown in Sect. 2, the claim for a strictly 
positive cosmological constant from SNIa data procceeds 
from an a priori homogeneity assumption. This, and the 
central role played in the whole cosmological field by the 
Cosmological Principle, is sufficient motivation to exam- 
ine to which extent this Principle can be observationally 
tested. 



As has been stressed by Ellis <| 1979[) . the situation 
is completly different within and outside our past light 
cone. For models exhibiting particle horizons, there are 
regions of the universe, outside our past light cone, from 
which we cannot receive any information at the present 
time. There is thus no way we can observationally verify 
the spatial homogeneity of such far out regions. 

The possibility of direct verification on our past light 
cone has been studied by Partovi and Mashhoon (j 1984) . 
These authors have explored the extent to which it is 
possible to distinguish between large scale homogeneous 
and inhomogeneous spherically symmetrical models, 
using magnitude-redshift data. Owing to the then state of 
observational technique, they concluded this distinction 
could not be performed. 

Analyses of galaxy redshift surveys have 
recently received increased attention. A con- 
troversy over whether the universe is smooth 
on large scale (|Guzzo 19971 
IMartinez et al. 19981 

Wu, Lahav fc Rees 1999[ ) or presents an 



Cappi et al 1998 



Scaramella ct al. 1998 



bounded fractal hierarchy ( Sylos Labini et al. 1 996, 
|Sylos Labini, Montuori fc Pietronero 1998 ) has devel- 
oped and waits for the next generation of wider and 
deeper galaxy catalogues which may provide a more 
conclusive answer (jMartinez 1999|) . 

The supernovae data have already been investigated 
for what they can tell about (in)homogeneity. Kim et 
al. (|1997p used the first seven SNIa discovered by the 
Supernova Cosmology Project at 0.35 < z < 0.65 and 
compared them to a nearby sample at z < 0.1 to declare 
the ruling out of the hypothesis of a locally underdense 
bubble. The nearby sample has also been examined for 
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evidence of a local "Hubble Bubble" by Zehavi et al. 
(1995J. The marginal signal identified by these authors 
is very sensitive to the cosmological model retained. 
Its value has been estimated for an Einstein-de Sitter 
universe, but would be less significant, and could be 
considered as included in the limits calculated by Kim et 
al., for the type of models proposed by R98 and P99. 

At another redshift extreme, Starkman et al. (1999) 
have shown that the segment of the universe sampled 
by the current supernovae data is not large enough 
to determine the overall properties of the expansion. 
But these authors only consider the Einstein-de Sitter 
homogeneous case to complete their caculations. 

It will however be shown, in Sect. 2, that in the near 
future, we would in principle, be able to probe homogene- 
ity on scales up to, at least z = 1. In principle here means 
provided the evolution effect comes under control, as well 
as any other bias or systematic effect. But, as long as the 
relevance of the homogeneous models used by R98 and 
P99 as a framework for their data analyses is not veri- 
fied, any claim to cosmological parameter measurements 
remains premature. Furthermore, a straight reading of the 
presently published results does not exclude a ruling out 
of the homogeneity hypothesis. It is moreover shown, in 
Sect. 3 and 4, that large scale inhomogeneity can mimic a 
cosmological constant in an homogeneous universe up to 
the precision achieved by current measurements. 

2. Magnitude-redshift relation and homogeneous 
models 

The luminosity distance Dl of a source is defined as the 
distance from which the radiating body, if motionless in 
an Euclidean space, would produce an energy flux equal 
to the one measured by the observer. It thus verifies 

l=— ^o". (1) 

L being the absolute luminosity, i.e. the luminosity in the 
rest frame of the source, and I the measured bolomet- 
ric flux, i.e. integrated over all frequencies by the observer. 

In Friedmann-Lemaitre-Robertson- Walker (FLRW) 
models, the distance measure at redshift z is a 
function of z and of the parameters of the model 
( Carroll, Press k Turner~l~992l) 



d l = s ( vT«T A 1 + z ') 2 ( 1 + 

-z'(2 + z')n A ]-^dz'\ (2) 



Hq being the current Hubble constant, Q,m, the mass den- 
sity parameter, and, £Ia being defined as 

where A is the cosmological constant, 
and with 

for £Im + > 1 i? = sin and k = 1 — Qm — 
for + f^A < 1 S = sinh and k = 1 — Qm — 
for VLm + = 1 S = I and k = 1 . 

The apparent bolometric magnitude m of a standard 
candle of absolute bolometric magnitude M, at a given 
redshift z, is thus also a function of z and of the param- 
eters of the model. Following Perlmutter et al. (Q~997), 
hereafter referred to as P97, it can be written, in units 
of megaparsecs, as 

to = M + 5logD L (z;fl M , Cl A ,H ) + 25 

= M + b\ogV L {z;n Mi n K ). (3) 

The magnitude "zero-point" M. = M — 5 log H + 25 
can be measured from the apparent magnitude 
and redshift of low-redshift examples of the stan- 
dard candles, without knowing Hq. Furthermore, 
V L (z; Q M , fi A ) = H D l (z;CIm,^a,Ho) depends on Q M 
and f^A with different functions of redshift. A priori 
assuming that FLRW models are valid to describe the 
observed universe, R98 and P99 thus proceed as follows 
to determine fijvf and 

A set of apparent magnitude and redshift measure- 
ments for low-redshift (0.004 z ^ 0.1) SNIa is used to 
calibrate Eq.([3]), and another set of such measurements 
for high-redshift (0.16 z ^ 0.97) is used to determine 
the best fit values of Qm and fiyjj. The results of these 
measurements are plotted in Fig. 4 and 5 of R98 and 
Fig. 1 and 2 of P99. The magnitude-redshift relation 
from the data and the theoretical curves obtained for 
different values of the cosmological parameters arc then 
compared and best-fit confidence regions in the parameter 
space arc plotted, see Fig. 6 and 7 of R98 and Fig. 7 of P99. 

It is convenient at this stage to make the following 
remarks. 

If one considers any cosmological model for which the 
luminosity distance is a function of the redshift z and 
of the other cosmological parameters of the model, and if 
this function is Taylor expandable near the observer, i.e. 

The methods used by the two survey teams, the stretch 
factor method for P99 and the multicolor light curve shape and 
template fitting methods for R98, differ in their conceptions 
and byproducts: P99 calibrate Ai and obtain direct fitting of 
£2m and J1a; R98 calibrate M, which yields an estimate for Ho- 
However, the following remarks apply to each of them. 
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around z = 0, the analysis of observational data at z < 1, 
in the framework of this model, can legitimately use the 
Taylor expansion 



D L (z;cp) = 
6 V dz 3 



z=0 



z=0 



1 

2+ 2 



dz 2 



+ 24 



1 /d 4 L>, 



dz 4 



z 

z=0 

4 i //V~5\ 



z 4 + 0(z 5 ), (4) 



z=0 



as, by definition of luminosity distances, Dl(z = 0) = 0. 
Here, cp denotes the set of cosmological parameters, 
pertaining to the given model, which can be cither 
constants, as in FLRW models, or functions of z, as in 
the example presented in Sect. 3 and 4. 

Luminosity distance measurements of sources at 
different redshifts z < 1 yield values for the different 
coefficients in the above expansion. Going to higher 
redshifts amounts to measuring the coefficients of higher 
power of z. For very low redshifts, the leading term is 
first order; for intermediate redshifts, second order. Then 
third order terms provide significant contributions. For 
redshifts approaching unity, higher order terms can no 
longer be neglected. 

Therefore, for cosmological models with very high 
(or infinite) number of free parameters, such data only 
provide constraints upon the values of the parameters 
near the observer. But for cosmological models with few 
parameters, giving independent contributions to each 
coefficient in the expansion, the method not only provides 
a way to evaluate the parameters, but, in most cases, to 
test the validity of the model itself. 



remarks, with D replacing T>, apply to the analysis and 
results of R98. 



In standard models, i.e. for f2 A = 



vf = ~(2-o M ) 



(9) 



If 



the 



analysis of the measurements gives 
V^' < 2 V L , it implies VLm < 0, which is physi- 
cally irrelevant. The model is thus ruled out. This is what 
happens with the SNIa data, and what induces R98 and 
P99 to postulate a strictly positive cosmological constant, 
to counteract the — £Im term. 

To test FLRW models with 11a ^ 0, one has to go 
at least to the third order to have a chance to obtain a 
result. This seems to be the order currently reached by 
the SNIa surveys. This last assertion rests on the two 
following remarks: 



1. As stressed in P97, R98 and P99, the well-constrained 
linear combination of fiji/ and Q\ obtained from 
the data is not parallel to any contour of constant 
current "deceleration'!! parameter qo = VLm/2 — ^a- 



As V {2) = V K l'{\ - q )/2, this implies that higher 
order terms effectively contribute. 

The contribution of the fourth and higher order terms 
is negligeable. One can easily verify that, for the 
higher redshifts reached by the surveys, the contribu- 
tion of the fourth order term does not overcome the 
measurement uncertainties, of the order of 5 to 10% 
(see Fig. 4 and 5 of R98 and Fig. 1 and 2 of P99). 



,(i) 



For FLRW models, precisely, one obtains from Eq.f2j) 

A ruling out of FLRW models with non zero cosmolog- 
(5) ical constant, at the third order level, i.e. due to a negative 
value for S]j}, would occur provided (see Eq.® to ([7])) 



L 2 



D 



(3) _ 



1 



M 4) EE - 



1 

24 



dD L 

dz 



c 



d 2 D L \ 



dz 2 
d 3 D, 



6 \ dz 3 



d^Dz 

dz 4 



z=0 



-mn 2 A 



12 ffiM £1 



M"A 



C 



4H t 

c 

8i?o 

^2 



(2-Q M + 2fi A ), 



(-2fi M -4ft A 



+^ + 40 A ), 
5c 



(0) 

4fi M fiA 

(7) 



72H 

Y2 



(8 A + 4n M fiA 



2 n 2 M 



80 A ) 



(8) 



1 - 



V 



(3) 



12) 



V 



(2)' 



< 0. 



(10) 



This cannot be excluded by the results of R98 and 
P99. It corresponds to the left part of the truncated 
best- fit confidence ellipsoidal regions in the 17m — 
plane of R98 Fig. 6 and 7 and P99 Fig. 7, and to the 
upper part of the error bars for the higher redshifts data 
in R98 Fig. 4 and 5 and P99 Fig. 1 and 2. 



These coefficients are independent functions of the 
three parameters of this class of models, Hq, Qm and 
With the method retained by P99 and described above, 
Hq is hidden in the magnitude "zero-point" M.. The 
coefficients T>^ of the expansion of T>l are thus functions 
of Qm and Q\ alone. However, the same following 



P99 propose an approximation of their results, which 
they write as 

0.8 M -0.6fi A » -0.2 ± 0.1, (11) 

2 A better name would be "acceleration" parameter, as qo < 
implies an accelerated expansion for a homogeneous universe. 
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which corresponds, in fact, to: 



V 



(2) 



V 



(1) 



4.2 



V 



(2) 



V 



(3) 



V 



(1) 



V 



(1) 



-1.8 ± 0.1. 



(12) 



It must be here stressed that the results published 
as primary by P99 and R98, under the form of best-fit 
confidence regions in the S]« — Qa plane, proceed from 
a Bayesian data analysis, for which a prior probability 
distribution accounting for the physically (in Friedmann 
cosmology) allowed part of parameter space is assumed. 
Results as given in the form of Eq. (fTTj) are thus distorted 
by an a priori homogeneity assumption, which would 
have to be discarded for the completion of the test here 
proposed. 

It may however occur that future, more accurate mea- 
surements yield values for the X^s verifying Eq. (|10[) with 
the > sign, which would correspond to a physically con- 
sistent positive r^Af • In this case, the FLRW models will 
have to be tested to the fourth order, i.e. with sources 
at redshifts nearer z = 1. A final test for the homogene- 
ity hypothesis on our past light cone would be a check of 
the necessary condition, obtained from Eq.© to ([8|), and 
which can be written as 



V 



(4) 



10 V 



(2)2 



20 V 



(2)3 



V 



(1) 



9 pCDa 



— VY> 



io vy'v 



(2W3) 



V 



(1) 



(13) 



The above described method only applies to data 
issued from z < 1 supernovae. If the ongoing surveys were 
to discover more distant sources, up to redshifts higher 
than unity, the Taylor expansion would no longer be valid. 

In practice, one will have to consider the Hubble 
diagram for the largest sample of accurately measured 
standard candles at every available scale of redshift from 
the lowest to the highest (e.g. Fig. 4 and 5 of R98 and 
Fig. 1 and 2 of P99). The FLRW theoretical diagram best 
fitting the data at intermediate redshift (z w 0.6 — 0.7 
for measurements with 5% accuracy, z w 0.8 for 10% 
accuracy) will be retained as candidate. If it corresponds 
to a negative value for SIm, the assumption of large scale 
homogeneity for the observed part of the universe will 
have to be discarded, and no higher redshift data will be 
needed, at least to deal with this issue. If it corresponds 
to a positive Qm, the diagram will have to be extended 
to higher redshift data. If these data confirm the best 
fit of this candidate model up to z w 0.8 — 0.9 (for 
measurements at some 5-10% accuracy), the homogeneity 
hypothesis (and the corresponding values of the model 
parameters) would receive a robust support (provided 



no unlikely fine tuning of another model of universe 
parameters). 

Another practical way of analysing the data is given 
by the very clever method described by Goobar and 
Perlmutter ( 1 99") j. hereafter refered to as GP. An appli- 
cation of this method to probe large scale homogeneity 
of the observed universe would, in principle, imply the 
measurement of the apparent magnitude and redshift 
of only three intermediate and high redshift standard 
candles. In practice, more would certainly be needed 
to smooth out observational uncertainties. The method 
applies to sources with any redshift values and runs as 
follows. 

Using Eq.([2]) and ([3]), one can predict the apparent 
magnitude of a source measured at a given redshift, in a 
peculiar FLRW model (i.e. for a given pair of values for 
VLm and J^a)- Fig. 1 of GP shows the contours of constant 
apparent magnitude on the f^A-versus-f^M plane for two 
sufficiently different redshifts. When an actual apparent 
measurement of a source at a given redshift is made, the 
candidate FLRW model selected is the one with values for 
VLm and J7a narrowed to a single contour line. Since one 
can assume some uncertainty in the measurements, the 
allowed ranges of JIm and £l\ are given by a strip between 
two contour lines. Two such measurements for sources at 
different redshift can define two strips that cross in a more 
narrowly constrained "allowed" region, shown as a dashed 
rhombus in Fig. 1 of GP. Now, if a third measurement 
of a source at a redshift sufficiently different from the 
two others is made, a third strip between two contour 
lines is selected. If this strip clearly crosses the previously 
drawn rhombus (and if the measurement uncertainties 
are kept within acceptable values), it provides support 
for the homogeneity hypothesis, with the above pointed 
out reserve. In case this strip clearly misses the rhombus, 
this hypothesis is ruled out, at least on our past light cone. 

If the latter happens to be the case, we shall need an 
alternative model to fit the data. An example of a model, 
able to fulfil this purpose without the help of a cosmolog- 
ical constant is proposed in the following section. 

3. Example: Lemaitre-Tolman-Bondi model with 
zero cosmological constant 

The complexity of the redshift-distance relation in 
inhomogeneous models and their deviation from the 
Friedmann relation have recently been stressed by Kurki- 
Suonio and Liang (fT992|) and Mustapha et al. (fT998 ). In 
a previous very interesting paper, Partovi and Mashhoon 
(198?|) have shown that the luminosity distance-redshift 
relation, in local models with radial inhomogeneities 
and the barotropic equation of state for the source of 
gravitational energy, cannot be distinguished from the 
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FLRW one, at least to second order in z. A special case 
of the solutions studied by these authors is used here to 
explore to which extent such inhomogeneous models can 
mimic homogeneity. 

A class of spatially spherically symmetrical solutions 
of Einstein's equations, with dust (pressureless ideal gas) 
as a source of gravitational energy, was first proposed by 
Lemaitre (|1933[) . It was later on discussed by Tolman 
(1934]) and Bondi ((1947), and became popular as the 
"Tolman-Bondi" model. In the following, it will be refered 
to as LTB model. 



2. with E(r) = 0, for all r 



R(r,t) = 



9GM(r) 



1/3 



3. with E(r) < 0, for all r 
GM(r) 



R 



t - t„(r) 



(19) 



-2E(r) 

GM(r) 
[-2E{r)} 3 / 2 



(1 — cosm), 



(20) 



Where the Cosmological Principle could 
to be ruled out, the LTB solution would ap- 
pear as a good tool for the study of the ob- 
served universe in the matter dominated region 
(jCelerier fc Schneider 19981 ISchneider fc Celerier 1999]) . 
It is used here as an example to show that a non van- 
ishing cosmological constant in a FLRW universe can 
be replaced by inhomogeneity with a zero cosmological 
constant to fit the SNIa data. 

The LTB line-element, in comoving coordinates 
(r, 9, ip) and proper time t, is 

ds 2 = -dt 2 + S 2 (r,t)dr 2 + R 2 (r,t)(d0 2 + sin 2 6df 2 ), (14) 
in units c = 1 . 

Einstein's equations, with A =0 and the stress-energy 
tensor of dust, imply the following constraints upon the 
metric coefficients: 

? ' 2 M) 



S\r,t) 



1 



R 2 (r, t) 



ATrp(r, t) 



R_ 

1 + 2E(r) 
GM(r) 
R(r,t) ~ 
M'(r) 



E(r), 



(15) 
(16) 
(17) 



R'(r, t)R 2 (r, t)' 

where a dot denotes differentiation with respect to t and 
a prime with respect to r, and p(r, t) is the energy density 
of the matter. E(r) and M(r) are arbitrary functions of r. 
E{r) can be interpreted as the total energy per unit mass 
and M(r) as the mass within the sphere of comoving 
radial coordinate r. 



One easily verifies that Eq. (|T6| possesses solutions for 
R{r, t), which differ owing to the sign of function E{r) and 
run as follows. 



1. with E(r) > 0, for all r 
GM(r) 



R 



t - t (r) 



2E(r) 
GM(r) 



[2£(r)] 3 / 



(coshu — 1), 
(sinh u — u). 



(18) 



where to(r) is another arbitrary function of r, usually 
interpreted, for cosmological use, as a Big-Bang sin- 
gularity surface H, and for which i?(r, t) = 0. One can 
choose to(r) = at the symmetry center (r = 0) by an 
appropriate translation of the t — const, surfaces and 
describe the universe by the t > to(r) part of the (r,t) 
plane, increasing t corresponding to going from the past 
to the future. 

The physical interpretation of R{r, t) must also be 
discussed. Bondi (1947J) presents it as the radial luminos- 
ity distance of a radiating source. But, with the definition 
used here (above Eq.{T])) and a corrected expression for 
apparent luminosity, it corresponds in fact to the luminos- 
ity distance divided by a factor (1 + z) 2 , as is shown below. 

Assuming that the light wavelength is much smaller 
than any reasonably defined radius of curvature for the 
universe (geometrical optics approximation), Kristian and 
Sachs (1965) established that the apparent intensity of a 
source, as measured at any point of one of its emitted light 
rays by an observer with proper velocity u a , is 



I = n(k a u a ) 2 



(21) 



where /i is a scalar, corresponding to the magnitude of 
the electromagnetic tensor components as measured by 
the observer. 



1 + 2 



From the definition of redshift, one obtains (|Ellis 1971[) 

{k a u a ) s 



(k a u a ) ' 



22) 



where k a is tangent to the null geodesies on which light 
travels, the subscripts s and o denoting values at the 
source and at the observer, respectively. 



In fact, the physical interpretation of this surface is im- 
proper. One can always consider that as the energy density 
increases while reaching its neighbourhood, radiation becomes 
the dominant energy component, pressure can no more be ne- 
glected, and the LTB model no longer holds. 
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For a measure realised at the source by an observer 
motionless in the rest frame of the source 



(*a«°)2 = 1, 

which gives 
H 



I = 



{1 + zf 



(23) 
(24) 



Assuming that the electromagnetic tensor for the light 
emitted by a distant source verifies Maxwell's equations 
for vacuum (jKristian fc Sachs 19661 Bondi HHZj), it fol- 
lows, for a radial measurement realised on any light ray 
by an observer located at the symmetry center (r = 0), 
that 



I = 



L 



4irR 2 (l + z) 4 



(25) 



The definition of luminosity distance given by Eq.([T]) 
yields the expression retained by Partovi and Mashhoon 
(1984), namely 



D L = {1 + zfR. 



(26) 



One thus sees that the luminosity distance Dl is 
a function of the redshift z and through R, of the 
parameters of the model: M(r), E(r) and to(r). 

A light ray issued from a radiating source with coordi- 
nates (t, r, 9, (p) and radially directed towards an observer 
located at the symmetry center of the model, satisfies, 
from Ea.([T4"|) and (fl"5j) . 



dt 



R'(r,t) 



zdr. 



(27) 



For a given function R(r,t), Eq.(|27|) possesses an infi- 
nite number of solutions of the form t(r), depending on the 
initial conditions at the source or on the final conditions 
at the observer. One can thus consider two rays emitted 
by the same source at slightly different times separated by 
r. The equation for the first ray can be written 

t = T(r). (28) 
The equation for the second ray is therefore 

t = T{r) + T {r). (29) 
Assuming: r(r) <C T(r) for all r, Eq. ([27f implies 

dT{r) _ R'[r,T(r)} 



dr v /l + 2£(r)' 

dr(r) _ r(r) 8Bf 

dr ~~ ~ y/l + 2E(r)~dT 



Mr)}, 



(30) 
(31) 



4 Bondi improperly writes: I — fi (U — in Bondi's nota- 
tion). 



which gives the equation of a ray fEq. (|3"0")) ) and the 
equation for the variation of r(r) along this ray (Eq. pTj) ). 

If one considers r(r) as the period of oscillation of some 
spectral line emitted by the source and r(0), its period as 
measured by the observer at r = 0, the definition of the 
redshift is 



r(r) = 



r(0) 
1 + z' 



(32) 



From this equation, another way of writing the rate of 
variation of r(r) along a ray is 



dr(r) 
dr 



•r(0) 



dz 



1 



dr (1 



(33) 



One can always choose z as a parameter along the null 
geodesies of the rays and obtain from the above equations 



dr 

dz 



dt 

dz 



Vl + 2E(r) 
(1 + z)R'[r,T(r)]' 



(34) 



Eq. ([27| therefore becomes 
R'[r,T(r)] 



(1 + z)R'[r,T(r)] 



(35) 



Eq. ([34|) and (f35|) form a system of two partial differ- 
ential equations of which each null geodesic is a solution 
starting with z at the source and finishing at the observer 
with z = 0. 

Successive partial derivatives of R with respect to 
r and t, and derivatives of E(r) with respect to 7', 
evaluated at the observer, contribute to the expression 
of the coefficients of the luminosity distance (Eq. (|2"6"|) ) 
expansion in powers of z. It is therefore interesting 
to note the behaviour of R(r,t) and E(r) near the 
symmetry center of the model, i.e. near the observer 
( [Humphreys, Maartens fc Matravers 1997| ): 



R(r,t) 



R'(0,t)r + O{r 2 ), 
E(r) = \E"{0)r 2 + 0(r 3 ). 



(36) 
(37) 



One thus easily sees that R, R, and higher order 
derivatives of R with respect to t alone, vanish at the 
observer, as do E and E 1 . 

Expressions for the coefficients of the luminosity dis- 
tance expansion naturally follow; after some calculations, 
one obtains 



D 



(i) 



R!_ 

R' : 



(38) 
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4 2) = 



1 R' ( R'R' R" R" 



2R> 



R ? R'R' R > 2 y ' 



(39) 



(3 ) _ 1 R! ( 1 R'R' , _ R' 2 R' _ R' 2 R' 



GR' 



R' 



+ 4 



R! 



R' 



R' 



+3- 



R"R> 
— — + 


R'R" 


R' 


■ 3 

R' 


R!" 


E"\ 


■ 2 R' 3 


7?' / 



R'R' 



(40) 



with implicit evaluation of the partial derivatives at the 
observer. 

Following Humphreys, Maartens and Matravers 
(1997), one can adopt a covariant definition for the Hub- 
ble and deceleration parameters of a spherically symmetric 
inhomogeneous universe. These authors (see also Partovi 
and Mashhoon (1984)) derive expressions for Ho and qo 
at the observer, in units c = 1: 



H 



1 



R' 



nW I Rf / ' 

i V / 

go = -H D {2) - 3. 



(41) 
(42) 



straightforward. 

As the mass function M{r) remains constant with 
time, it can be used to define a radial coordinate r: 
M(r) = Mqt 3 , where Mq is a constant. 

With the covariant definition for Hq above mentioned 
(Eq.gT])), the V^s, as derived from Eq.® to (gDJ), can 



thus be written, in units c ; 

D w - _L 



1, in the form 



D 



(2) _ 



1 



D 



(3) 



4ff 



1-6- 



^ 9GAIq * 



(43) 



(44) 



/ 3 



^(0) 



^ 2 (0) 



^ 9GM ' 



^ 9GM \ 



/ 3 



(45) 



with the previously indicated choice to(0) = 0, where t p is 
the time-like coordinate at the observer. It is convenient 
to note that t p is not a free parameter of the model, 
since its value proceeds from the currently measured 
temperature at 2.73 K (Celerier & Schneider 1998). 



Substituting into Eq. (|3"8"|) to (|4H|) . it is easy to see 
that the expressions for the FLRW coefficients in the 
expansion of £>£ in powers of z can mimic LTB with 
A = ones, at least to third order. This is straightfor- 
ward for D [ l ] . The case of Z?£ 2 ' is discussed at length in 
Partovi and Mashhoon (1984). For higher order terms, 
it implies constraints on the LTB parameters, which 
will be illustrated below on the peculiar example of flat 
models. In fact, owing to the appearance of higher order 
derivatives of the parameter functions in each higher 
order coefficient, LTB models are completly degenerate 
with respect to any magnitude-redshift relation, while 
FLRW ones, of which the parameters are constants, are 
more rapidly constrained and cannot thus fit any given 
relation, when tested at sufficiently high redshifts. 



4. Illustration: flat LTB (A = 0) models 

To illustrate which kind of constraints can be imposed 
on LTB parameters by current observational results, the 
peculiar case of spatially flat LTB (A = 0) models is 
analysed here. 

Spatial flatness is a property of the subclass of LTB 
models verifying E(r) = (Bondi 1947). In this case, 
the expression for R is given above by Eq. (fT9)) and the 
calculation of the successive derivatives of i?, contributing 
to the expressions for the expansion coefficients, is 



A comparison with the corresponding FLRW coeffi- 
cients gives the following relations: 



M 



1 



9 
"4 

1 

' 2 

9 



15 

y 



t'o 2 (0) 



3 f 3 
Op 



(^) 
t'oiO) 

/ 9GM \ 1 



3 f 3 



15 

T 



' 9GMn ^ 
, 2 / 



^ 2 (0) 



(46) 



^ 9GM ' 



(47) 



The above Eq. (|47[) implies that a non vanishing 
cosmological constant in a FLRW interpretation of data 
at z < 1 corresponds to a mere constraint on the model 
parameters in a flat LTB (A = 0) interpretation. 

Any magnitude-redshift relation, established up to the 
redshifts and with the precisions achieved by current mea- 
surements, i.e. at third order level, can thus be interpreted 
in either model. For instance, the latest results published 
in P99, and given under the form of Eg. pip , correspond, 
in a flat LTB (A = 0) interpretation, to 



4.3- 



/ 9GAf -j 3 
-1 ± 0.1. 



3.75- 



^ 2 (0) 



1.125- 



i'o'(Q) 



' 9GM * 
v 2 / 



(48) 
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Such a result would imply a negative value for at least 
one of the two quantities t' (0) or £q(0), which would be an 
interesting constraint on the "Big-Bang" function in the 
observer's neighbourhood. For instance, a function to(r) 
decreasing near the observer would imply, for a source at 
a given z < 1 , an elapsed time from the initial singularity 
that is longer in an LTB model than in the corresponding 
FLRW one, i.e. an "older" universe 0. A decreasing to(r) 
has thus an analogous effect in a LTB universe to a positive 
cosmological constant in a FLRW one. They both make 
the observed universe look "older" . 

5. Conclusion and discussion 

The interpretation of recently published data from high 
redshift SNIa surveys, devoted to the measurement of the 
magnitude-redshift relation, has been re-examined with 
no a priori idea about which model would best describe 
our universe, at least up to redshifts z ~ 1. 

It has been shown that a straight reading of these 
data does not exclude the possibility of ruling out the 
Cosmological Principle. 

A method of testing large scale homogeneity on our 
past light cone, provided sufficiently accurate data from 
"standard candles" at redshifts approaching z = 1 would 
be available, has been proposed. This method could be 
applied using the SNIa data, provided every source of 
potential bias or systematic uncertainties would be cor- 
rectly taken into account. An evolution of the progenitors 
is generally considered as one of the more likely sources 
of potential systematic errors in the analysis of the super- 
novae data. Mustapha et al. (1997) have made the point 
that, given isotropic observations about us, homogeneity 
cannot be proved without either a model independent 
theory of source evolution, or distance measures that are 
independent of source evolution. If a perceptible evolution 
effect was to be put forward, its impact on the SNIa data 
would thus have to be evaluated with a model indepen- 
dent method, before using them as a test for homogeneity. 

An example of an inhomogeneous model with a zero 
cosmological constant reproducing the current observa- 
tions has been given. However, a vanishing cosmological 
constant is a minimal feature here imposed on the 
model only to prove that A is not necessary to fulfill the 
observational constraints. A non zero A would only add 
a new free parameter in the LTB equations and would 
not alter the following primary statement: the presently 
published SNIa data can be interpreted as implying 
cither a strictly positive cosmological constant in an 
homogeneous universe or large scale inhomogeneity with 
no constraint on A. 



The choice of over-simple inhomogeneous models 
as examples can, of course, be questioned. The most 
naive assumptions of spherical symmetry and centered 
observer, which have been discussed in other works (see 
e.g. ICelerier fe Schneider 1998[ ICelerier 1999[) . have been 
mainly retained for simplification purpose, but do not 
distort the above conclusions. It could be actually inter- 
esting to probe the possibility of an off-center observer 



In fact, see footnote 3, a lower temperature. 



(see e.g. [Humphreys, Maartens fc Matravers 1997[ ) with 
the SNIa data. Spherical symmetry is grounded on the 
observed quasi-anisotropy of the CMBR temperature, 
and can thus be considered as a sufficiently good working 
approximation. 

The possibility of testing finer features of the model 
depends on observational improvements. But before 
aiming at such ambitious goals, it is of the utmost 
importance to complete the large scale homogeneity test. 
A previous increase in the number and measurement 
accuracy of the candidate "standard candles" at very 
high redshifts is therefore urgently needed, for progress 
in both fundamental issues of the Cosmological Principle 
and of the cosmological constant. 

It has been shown however that LTB models are 
highly degenerate with respect to any magnitude-redshift 
relation, but this is not the case for FLRW models. 
The best way to prove large scale inhomogeneity would 
therefore be to disprove homogeneity. Conversely, the 
best way to prove a non zero cosmological constant 
would be to prove large scale homogeneity. But, as the 
Friedmann distance-redshift relation is a necessary but 
not a sufficient condition for homogeneity, its observed 
verification would not, in principle, be enough to support 
the Cosmological Principle. Even if this would imply 
a fine tuning of its parameters, the possibility for an 
inhomogeneous universe to mimic such a relation could 
not be excluded. 



To consolidate the robustness of the future magnitude- 
redshift tests, it would therefore be worth confronting 
their results with the full range of available cosmological 
data. The most recent attempt to do so is the work by 
Bahcall et al. |1999) who conclude that the analysed 
observations can be considered as consistent with the 
standard Big-Bang picture of the expansion of the 
universe, but do not discuss the homogeneity hypothesis. 
A cross-check with a model independent analysis of the 
CMBR anisotropy would, for instance, likely yield a sub- 
stantial improvement, as the CMBR data are known to 
provide orthogonal constraints on the models parameters. 
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